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We show how growth by agglomeration can be described by means of algebraic or differ- 
ential equations which determine the evolution of probabilities of various local configu- 
rations. The minimal fluctuation condition is used to define vitrification. Our methods 
have been successfully used for the description of glass formation. 



1. Introduction 

In a series of papers published during the past ten years 

ft S, I B) 

, new models 

of growth by agglomeration of smaller units have been elaborated, and applied to 



many important physical systems, S.U£h quasicrystals (El), fullerenes (Q, 0), and 
oxide and chalcogenide glasses, (□, u,tL2l, til). Here we shall present the main ideas 
on which these models are based, and briefly discuss the latest developments. 

In order to make our presentation concise, the example we choose is the sim- 
plest covalent network glass known to physicists, the binary chalcogenide glass 
As x Se(\- x ), where x is the concentration of arsenic atoms in the basic glass-former, 
which in this case is pure selenium. The generalization to other covalent networks, 
e.g. Ge x Se(\^ x ), is quite straightforward. These glasses (in the form of thin and 
elastic foils) are used in photocopying devices. 

Whether the formation of a solid network of atoms or molecules occurs in a 
more or less rapidly cooled liquid, or as vapor condensation on a cold support, the 
most important common feature of these processes is progressive agglomeration of 
small and mobile units (which may be just single atoms, or stable molecules, or even 
small clusters already present in the liquid state) into an infinite stable network, 
whose topology can no longer be modified unless the temperature is raised again, 
leading to the inverse (melting or evaporation) process. 

To describe such an agglomeration with all geometrical and physical parameters, 
such as bond angles and lengths, and the corresponding chemical and mechanical 



*Email: rk@ccr.jussieu.fr 



1 



6, 2008 18:9 WSPC/Guidclines tianjin3 



2 RICHARD KERNER 

energies stored in each newly formed bond, is beyond the possibilities of any rea- 
sonable model. This is why stochastic theory is an ideal tool for the description 
of random agglomeration and growth processes. Instead of reconstructing all local 
configurations, it takes into account only the probabilities of them being found in 
the network, and then the probabilities of higher order, corresponding to local cor- 
relations. This is achieved by using the stochastic matrix technique. A stochastic 
matrix M represents an operator transforming given finite distribution of proba- 
bilities , [pi,p2, ...,pn] ! into another distribution of probabilities, [pi,f>2) ■■■iPn]- ^ 
follows immediately that such a matrix must have only real non-negative entries, 
each column summing up to 1. 

The algebraic properties of such matrices are very well known. The main fea- 
ture that we shall use here is the fact that any stochastic matrix has at least one 
eigenvalue equal to 1. The remaining eigenvalues have their absolute values always 
less than 1. This means that if we continue to apply a stochastic matrix to any 
initial probability distribution, after some time only the distribition corresponding 
to the unit eigenvalue will remain, all other contributions shrinking exponentially. 
This enables us to find the asymptotic probability distribution. 

In what follows, we identify these probability distributions with stable or meta- 
stable states of the system, fixing the statistics of characteristic sites in the network. 
Taking into account Boltzmann factors (with chemical potentials responsible for the 
formation of bonds), we are able to find the glass transition temperature in various 
compounds. In particular, oneJs able to predict the initial slope of the curve T g {c), 

1. e. the value of {dT g /dc) c=0 (0, 0). 

2. Stochastic matrix describing cluster agglomeration 

Consider a binary selenium-arsenic glass, in which selenium is the basis glass former, 
and arsenic is added as modifier (although its concentration can be as high as 
30%). The chemical formula denoting this compound is As c Se(\- c ), where c is the 
As concentration. In a hot liquid, prior to solidification, the basic building blocks 
that agglomerate are just selenium and arsenic atoms, indicated respectively by 



appear everywhere, growing by agglomeration of new atoms on their rim. Consider 
a growing cluster: one can distinguish three types of situations (we shall call them 
"sites") on the cluster's rim. The concentration of free As atoms in the liquid will 
be called c and that of Se, (1 — c). 

We should stress here that the mathematical model we propose to analyze is 
quite far from reality in the particular case of As — Se binary glass, as has been 
shown in (U), because some of the AS atoms are five-coordinated. Our model 
gives better predictions for the Ge — Se binary glasses. Our aim here is to expose 
the basic theoretical concepts rather than get very precise predictions; for this, 
one needs more sophisticated models, taking into account various possibilities of 
different agglomeration modes, e.g. ring-forming and coordination changes. 




When the temperature goes down, clusters of atoms start to 
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Two choices are possible for constructing the states and transition matrix (see 
(lj)). There are three possible kinds of sites: a selenium atom with one unsaturated 
bond, and an As atom presenting one or two free bonds; these are indicated by x 
= o — , y = • and z = • — . To each site one of the two basic cells can attach itself, 
reproducing one of the initial configurations, in the specific combinations shown in 
the next column of the Figure 1. The attachment of one single basic cell, or the 
saturation of one single bond, is a step in the evolution. In the second choice, each 
step is obtained by the complete saturation of all the bonds at the rim, so that only 
two types of sites (denoted by x and y) are seen on cluster's rim, assum ing that 
the growth is of dendritic type (no small rings present). It can be shown (EJ) that 
the two approaches lead to the same results, which may be considered as a proof 
of the ergodicity of the proposed model. We shall choose the second version of the 
model for the sake of simplicity. In this case, we can take into account only the x 
and y-type sites, because the z-type sites transform after the next agglomeration 
step into an x or y type site. The elementary step in the agglomeration process, 
described by he transition matrix, corresponds now to the complete saturation of 
all the available free bonds on the rim. This is represented in Figure 1 : 

X 2 (1 - c) e- e 

y See"' 



-A 




2 X 4 (1 - cf e- 2 " 

— <^~~ x + y 12 c (1 - c) e""- a 

2 y 9 c 2 e' 2a 

Figure 1: States, steps and un-normalized probability factors . 

Observing that from the site z only the sites of x and y type can be produced, 
we can forget it and consider the dendritic growth with only two types of sites 
appearing all the time. Given an arbitrary initial state (p x ,P y ), the new state results 
from taking into account all possible ways of saturating the bonds of the previous 
state's sites by the available external atoms. The un-normalized probability factors 
arc displayed in the Figure. The non-normalized probability factors can be arranged 
in a matrix 



2(l-c)e- £ 4(l-c) 2 e" 2,) 
8(1 - c) 2 e- 2r i + 12c(l - c)e- r i- a 12c(l - c)e~''- a + I8c 2 e~ 2a 

The normalized transition matrix is written as 



(1) 



M= , M XX M X y \ / M XX 1 - Myy , 

\My X Myy) \ 1 - M xx Myy ' 
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where the entries are obtained by normalizing the columns of the matrix (Q). 

M « = 2(1-^+30 ' and M ^= 2(l- 3 T + 3c, < 3 > 

where we have introduced the abbreviated notation £ = e 1 ' -6 and /i = e 1 ^". 

The eigenvalues of this matrix are 1 and M xx — M yy = M xy — Myx. and the 
stationary eigenvector is 



P?\ 1 (M xy 

P™J M xy + M yx \M yx 



(4) 



It can be seen from Figure 1 that on the surface of an average cluster, p x is the Se 
concentration and p y is the As concentration. Now, the high homogeneity exhibited 
by known glass structures suggests that even in relatively small clusters, deviations 
from the average modifier concentration c must be negligible. Thus, in the bulk, 
the As concentration should be equal to c. Therefore, the condition of minimal 
fluctuations in the bulk concentration can be interpreted as the glass transition 
condition. This means that the asymptotic state is fixed by the external concen- 
tration, therefore the above eigenvector must be equal to the average distribution 
vector (1 — c, c). The solutions are c = 0, c = 1 and the nontrivial one 

M yx 6 - 4£ 



M xy + M yx 12 - 4f - 9/j 

This equation can be checked against experiment. For example, we can evaluate 
the derivative = (^f ) for a given value of c. In particular, as c — -> 0, we can 
neglect the As- As bond creation (equivalent to putting \i — in (||)), to get 



dT 
~dc~ 



T 

± go 



ln(3/2) 



(where T go is the glass transition temperature of pure Se). This is the present-case 
expression of the general formula given by the stochastic approach, the fraction 
(3/2) being replaced by (m'/m), where m and m' are the valences of the basic glass 
former and o f th e modifier) , remaining in very good agreement with the experimen- 



tal data (see M, M, ED) 



3. Low concentration limit. 

The above scheme can be easily generalized to the case of arbitrary valence, say 
rriA and tub- In that case, the stochastic 2x2 matrix has the same form as (^J), 
but with the entries given by 

7714(1 — c)£ irijifl-c) 

M XX = 1 ~M yX = M / C , M Xy = 1 — Myy = A \ ' , 

777,^(1 — c)£ + rase niA(l — c) + rriBCfJ, 
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The asymptotic probability has the same form as before, as well as the zero fluctu- 
ation condition relating c with T (interpreted as the glass transition temperature) . 
The derivative of c with respect to the temperature T gives the "magic formula" 

dT T [(1 - + (1 - ^ M )P U 

where we used the fact that ^ = — A £ ln£ , and -J£ = \xln[L. This defines 
the slope of the function T ff (c), which is an important measurable quantity : 

dT g [(l-^) + (l-^M)] 2 

The initial slope, at c = 0, is of particular interest. Its expression is very simple, 
taking into account that when c = 0, we have also £ = ^J, which leads to 



dT g 



dc 



V ra a J 



Its value has been checked against the experiment very successfully, in more than 30 
different compounds. In some cases the formula does not seem to work well; usually 
it comes from the change of valence of certain atoms provoked by the influence of 
the surrounding substrate. 

One could be worried about the apparent singularity in this formula when tua = 
tob, i-e. when one deals with a mixture of two different glass formers with the same 
coordination number. It is not difficult to show that also in such c\ ('MSG cL r easonable 
limit can be defined, as has been recently suggested by M.Micoulaut (113). As a 
matter of fact, suppose that the glass transition temperature of the pure glass- 
former A is T g o, and that of the pure glass- former B is T g \. We can re- write our 
minimal fluctuation condition in a very symmetric manner, invariant with respect 
to the simultaneous substitution rriA ^ tub, c <-» (1 — c) and £ <-» fi : 

c (l- c )[(l- c )(l-™±0-c(l-— M)]=0 (9) 
m B m A 

Obviously, the "pure states" c = or c = 1 represent stationary solutions of 
and can be factorized out. The non-trivial condition for the glass forming is thus 

(l-c)[l-^]-c[l-^ M ]=0 (10) 
m B m A 

Now, using the limit conditions at c — > 0, T g = T g o and c — > 1, T g = T g ±, 
and introducing the generalized Boltzmann factors with the energy barriers for 
corresponding bond creations as Eaa, Eab and Ebb, we can write 

E A b ~ E A a — kT g0 ln( — — ) , Eab - Ebb = kT gl ln( — — ), (11) 
m A m B 

so that the expressions £ and \i at the arbitrary temperature T can be written as 



1 O 



,ms. r 90 VAB-^BB ,^gl mA- 



£(T) = e =(_S.)T- ; /i(T) = e ^ - = ( — )~ ■ (12) 

m A m B 
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Substituting these expressions into (]?]) and taking the limit c — > 0, we get 

— (13) 



dT 9| _T fl0 [l-(^)" 



It is easy to see now that even when tua = tub, this formula has a well defined 
limit. Indeed, if we first set = 1 + e, and then develop the numerator and the 
denominator of the above equation in powers of e, then in the limit when e — > we 
arrive at a simple linear dependence which is in agreement with common sense and 
with experiment as well, namely 

^\ c =o=T gl -T g0 (14) 

This formula is also confirmed by many experiments, e.g. performed on selenium- 
sulfur mixtures (where ttia = tiib = 2). The deviations from the linear law (|l4j) 
observed in the Se — Te binary glass are explained by the fact of the chemical 
properties of tellurium, which changes its valence from 2 to 3 in presence of selenium. 

4. The effect of rapid cooling 

An interesting extension of this model is obtained when we take into account the 
effects of rapid cooling, i.e. when the time derivative of the temperature can no 
longer be neglected. The treatment of this problem was suggested in (E3), and has 
been solved quite recently (EJ). 

Consider the agglomeration process defined by the above stochastic matrix, 
p' = Alp, with p representing a normalized column (a "vector") with two entries, 
p x and p y = 1 — p x . After one agglomeration step, representing on the average one 
new layer formed on the rim of a cluster, we can write 

Ap =p' —p = (M — l)p (15) 

Let us introduce a symbolic variable s defining the progress of the agglomeration 
process; obviously, s(t) should be a monotonically increasing function during the 
glass transition. If the temperature variation is so slow that the derivative dT/dt = 
(dT/ds)(ds/dt) can be neglected (which is often called the annealing of glass), the 
master equation of our model can be written as 

Ap=^As = (AI-l)pAs 
os 

where the variation As represents one complete agglomeration step. If we want to 
use real time t as an independent parameter, we should write 

f=^ = r~ 1 ^=-(M-l)p (16) 
dt Asdt As t v ' y ' 

We have introduced here the new entity r = (ds/dt)^ 1 which can be interpreted as 
the average time needed to complete a new layer in any cluster, or alternatively, the 
time needed for an average bond creation. Now, if the temperature varies rapidly 
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enough, the matrix M can no longer be considered as constant. The equation ( |T6| ) 
must be modified according to the well known "moving target" principle. That is, 
the total derivative of p with respect to t should read: 



d P I** i\ ds _ dM dT _ dp 
_ = (M -l)-p + — -„ = - 



(M - 1) + q 



dM 
~df~ 



P 



(17) 



where we supposed linear dependence of the temperature on time, so that the 
derivative dT/dt can be denoted by constant cooling rate q. In the two-dimensional 
case only one component of p is independent, because p x +p y = 1. Let us choose p y 
(whose asymptotic value should be equal to c) as independent variable. Then (|l7j ) 
will reduce to the single equation : 



dpy 
dt 



1 



(Myy ~ l) Py + My X (l - P y) 



q 



8M„ 



yy 



dT 



Py 



9M yx 
dT [ 



Py) 



(18) 



where we have used the fact that p x — 1 — p y , M xx = 1 — M yx and M yy = 1 — M xy . 

What remains is just simple algebra. After a few operations we find the asymp- 
totic value of p y , 



denoted obtained when we set dp y /dt = 
M vx + rq( 

Pv 



'V 



{M xy + Myx) + rq{ 



dT 



d(M*y+M y3C ) 
dT 



(19) 



As in the former case, we define the glass transition temperature by solving the zero- 
fluctuation condition p^ = c. The quasi-equilibrium condition thus obtained can 
be written in a form displaying an apparent symmetry between the two ingredients 
("A" and "B") of binary glass. As in the previous case (when q = 0), the limit 
values c = and c = 1 represent stationary solutions, which is obvious (no local 
fluctuations of concentration c are possible when there is no ingredient other than 
A or B atoms alone). After factorizing out c(l — c), we get 



m B 



m A 



7TIa(1 — c)£ + tubc ?tm(1 — c) + tubc/i 



rq 



■ m A m B 



c^tln/j, 



(l-c)£/n£ 



where we have used the fact that d ^ L ^' > 



[m^(l — c) + m B c^] 2 [toa(1 - c) £ + mgc] 2 



(20) 



The 



d(ln fj.) _ In fi 

T ' dT ~ T 

above formula seems quite cumbersome, but it become much simpler in the low 

Close to c 



concentration limit, c 



m B 
mA 



we get 
r q tub 



ln( = Q 



(21) 



T m,A 

1 one gets the same formula switching tua with 



(quite obviously, in the limit c 
m B and replacing £ by /i). Replacing £ by the expression (12), we arrive at : 



( mB 
\m A 



0. 



(22) 
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It is easy to see that independently of the ratio m^/m^, for temperatures T above 
T g o we must have q < 0, and vice-versa, during rapid cooling the glass transition 
occurs at the temperature T > T g0 . 

The dimensionless combination (r q) /T defines the quenching rate as the prod- 
uct of (l/T)(dT/dt) = d(lnT)/dt by the time constant r, characterizing the kinetics 
of the agglomeration process, i.e. the average time it takes to create a new bond. It 
may depend weakly on the temperature, but for the sake of simplicity suppose it 
is constant. It can be determined by comparing formula (p2|) with the experimen- 
tal data. To take an example, let us again consider the selenium-arsenic glass at 
c — > (almost pure selenium with a small addition of As). We know that in this 
case T g — > T g o — 318°K. The formula (|2^) then gives the quasi-linear dependence 
of AT = T - T gQ on the quenching rate q: for T g = 328°K (i.e. AT = 10°i^) we 
get rq = -10.38; for T g = 338°^ (i.e. AT = 20°K) we get rq = -21.51; for 
T g = M8°K (i.e. AT = 30°^) we get rq = -32.26, and so forth. 

Finally, if we want to establish the formula for a pure glass-former, without any 
modifier, we should take the limit (rUA/ms) — > 1 and /i — > £; we then get 

^^ + (^)^ = or T-T = AT g = -(rq)^f. (23) 

Eventually, the deviations from this simple dependence may indicate that the char- 
acteristic time t depends on T. This can shed more light on the agglomeration 
kinetics in various glass-forming liquids. More details can be found in (El), (§3) 

Acknowledgements 

Enlightening discussions with R.Aldrovandi, R. A. Barrio and M.Micoulaut are grate- 
fully acknowledged. 

References 

1. R. Kerner, Physica B215 (1995) 267 

2. R.Kerner and M.Micoulaut, Journ. of Physics: Cond. Matter, 9, 2551-2562, (1997). 

3. R.Kerner, A theory of glass formation, in Atomic diffusion in amorphous solids, M. 
Balkanski and R.J. Elliott, eds., World Scientific, p.25-80 (1998) 

4. R.Kerner, The principle of self-similarity and its applications to the description of non- 
crystalline matter, Proceedings of the Workshop in Cocoyoc, 1997, J.L. Moran-Lopez 
ed., Plenum Press, 323-337, (1998). 

5. D.M. dos Santos-Loff and R.Kerner, Journal de Physique I (4), 1491-1511, (1994 

6. R. Kerner, K.H. Bennemann, K. Penson, Europhysics Lett. 19 (5), 363-368 (1992). 

7. R. Kerner, Computational Materials Science, 2 p. 500-508, (1994); 

8. D.M. dos-Santos-Loff, R.Kerner and M.Micoulaut, Journ. of Phys. C, 7, 8035-8052, 
(1995). 

9. R.A.Barrio, J.P.Duruisseau and R.Kerner, Phil. Magazine B, 72 (5), 535-550, (1995). 

10. R.A.Barrio, R.Kerner, M.Micoulaut, G.G.Naumis, Journ. of.Phys.: Condensed Mat- 
ter, 9 p.9219-9234 (1997) 

11. R.Kerner, G.G.Naumis, Journal of Non-Cryst. Solids 231 p.111-119 (1998) and 
R,Kerner, G.G.Naumis, Journ. of Physics: Cond. Matter 12 (8), p. 1641-1648, (2000). 



6, 2008 18:9 WSPC/Guidelines tianjin3 



Stochastic agglomeration and growth 9 

12. R.Kerner, Journ. of Non-Cryst. Solids, 182, p.9-21 (1995) 

R.Kerner, M.Micoulaut, Journ. of Non-Cryst. Solids, 210, p. 298-305 (1997) 

14. R.Aldrovandi and R.Kerner, in "New Symmetries and Integrable Models", Proceedings 
of XIV Max Born Symposium, eds. A. Frydryszak, J. Lukierski and Z. Popowicz, World 
Scientific, p. 153-169, (2000) 

15. PBoolchand and W.J.Bresser, Phil.Mag.B 80, 1757-1772 (2000) 

16. Boolchand P., Feng X. and Bresser W.J., Rigidity transition in binary Ge — Se glasses, 
Journ. Non-Cryst. Solids 293, p.348 (2001) 

17. D.G.Georgiev, PBoolchand and M.Micoulaut, Phys.Rev.B 62, R9228 (2000) 

18. R.Kerner, Mathematical models of glass formation , Proceedings of the conference 
"Glasses and Solid Electrolytes " (St. Petersburg, May 1999), Glass Physics and Chem- 
istry, 26 (4), p.313 - 324 (2000) 

19. M.Micoulaut, private communication (2001) 

20. R.Kerner, M.Micoulaut, C. R.Acad. Sci.Paris, t. 315, Ser.II, p. 1307-1313 (1992) 

21. R.Kerner and J.C.Phillips, Solid State Communications, 117, p.47-56 (2000); 
R.Kerner, in preparation. 



